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Abstract 

We define a new diffusive matrix model converging towards the /3-Dyson Brownian 
motion for all f3 e [0,2] that provides an explicit construction of /3-ensembles of random 
matrices that is invariant under the orthogonal/unitary group. We also describe the 
eigenvector dynamics of the limiting matrix process; we show that when /? < 1 and that 
two eigenvalues collide, the eigenvectors of these two colliding eigenvalues fluctuate 
very fast and take the uniform measure on the orthocomplement of the eigenvectors 
of the remaining eigenvalues. 
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1 Introduction 

It is well known that the law of the eigenvalues of the classical Gaussian matrix ensembles 
are given by a Gibbs measure of a Coulomb gas interaction with inverse temperature (3 = 1 
(resp. 2, resp. 4) in the symmetric (resp. Hermitian, resp. symplectic) cases; 



Such measures are associated with symmetric Langevin dynamics, the so-called Dyson 
Brownian motion, which describe the random motion of the eigenvalues of a symmetric 
(resp. Hermitian, resp. symplectic) Brownian motion. They are given by the stochastic 
differential system 



with iid Brownian motions These laws and dynamics have been intensively studied, 
and both local and global behaviours of these eigenvalues have been analyzed precisely, 
starting from the reference book of Mehta [9]. 

More recently, the generalization of these distributions and dynamics to all /? > 0, the 
so-called /3-ensembles, was considered. As for f3 = 1,2,4, the Langevin dynamics converge 
to their unique invariant Gibbs measure Pp as times goes to infinity. Indeed, the stochastic 
differential system under study is a set of Brownian motions in interaction according to a 
strictly convex potential. Thus, one can then show by a standard coupling argument that 
two solutions driven by the same Brownian motion but with different initial data will soon 
be very close to each others. This entails the uniqueness of the invariant measure as well 
as the convergence to this Gibbs measure. It turns out that the case f3 £ [0, 1) and the 
case (3 £ [1, oo) are quite different, as in the first case the eigenvalues process can cross 
whereas in the second the repulsion is strong enough so that the eigenvalues do not collide 
with probability one in finite time. However, the diffusion was shown to be well defined, 
even for f3 < 1, by Cepa and Lepingle [4], at list once reordered. 

The goal of this article is to provide a natural interpretation of /3-ensembles in terms of 
random matrices for (3 6 [0,2]. Dumitriu and Edelman [6] already proposed a tridiagonal 
matrix with eigenvalues distributed according to the /3-ensembles. However, this tridiag- 
onal matrix lacks the invariant property of the classical ensembles. Our construction has 
this property and moreover is constructive as it is based on a dynamical scheme. It was 
proposed by JP Bouchaud, and this article provides rigorous proofs of the results stated 
in pp. The idea is to interpolate between the Dyson Brownian motion and the standard 
Brownian motion by throwing a coin at every infinitesimal time step to decide whether 
our matrix will evolve according to a Hermitian Brownian motion (with probability p) 
or will keep the same eigenvectors but has eigenvalues diffusing according to a Brownian 
motion. When the size of the infinitesimal time steps goes to zero, we will prove that the 
dynamics of the eigenvalues of this matrix valued process converges towards the /3-Dyson 
Brownian motion with (3 = 2p. The same construction with a symmetric Brownian motion 
leads to the same limit with {3 = p. This result is more precisely stated in Theorem 12.21 
We shall not consider the extension to the symplectic Brownian motion in this paper, but 



dp ^ x ) = -yp n i a * - e a? • 



i<3 
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it is clear that the same result holds with f3 = 4p. Our construction can be extended to 
other matrix models such as Wishart matrices, Circular and Ginibre Gaussian Ensembles 
and will lead to similar results. 

We thus deduce from our construction that /3-ensembles can be interpreted as an in- 
terpolation between free convolution (obtained by adding a Hermitian Brownian motion) 
and standard convolution (arising when the eigenvalues evolve following standard Brow- 
nian motions). It is natural to wonder whether a notion of /3-convolution could be more 
generally defined. 

Moreover we shall study the eigenvectors of our matrix- valued process. In the case 
where j3 > 1, their dynamics is well known and is similar to the dynamics of the eigen- 
vectors of the Hermitian or Symplectic Brownian motions, see e.g. [2]. When /3 < 1 the 
question is to determine what happens at a collision. It turns out that when we approach 
a collision, the eigenvectors of the non-colliding eigenvalues converge to some orthogonal 
family B of d — 2 vectors whereas the eigenvectors of the colliding eigenvalues oscillate 
very fast and take the uniform distribution on the ortho-complement of B, see Proposition 



2 Statement of the results 

Let Ti^ be the space of d x d symmetric (respectively Hermitian) matrices if j3 = 1 (resp. 
/3 = 2) and be the space of d x d orthogonal (respectively unitary) matrices if f3 = 1 
(resp. = 2). 

We consider the matrix-valued process defined as follows. Let 7 be a positive real 
number and Mq £ with distinct eigenvalues Ai < A2 < • • • < A^. For each n G N, we 
let (e£)fc g N be a sequence of i.i.d {0, l}-valued Bernoulli variables with mean p in the sense 
that 

P[e£ = l]=p=l-P[e£ = 0]. 

Furthermore, for t ^ 0, we set e™ := c? nt y 

In the following, the process (t))t ^ will denote a symmetric Brownian motion, 
i.e. a process with values in the set of d x d symmetric matrices (respectively Hermitian if 
P = 2) with entries H^-{t),t ^ 0, i ^ j constructed via independent real valued Brownian 
motions (Bij , By , 1 ^ i ^ j ^ d) by 

H^(t\ = I B ^ + W ~ VBv® ' lii< i (OK 
« W \ y/2Ba(t) otherwise 1 ' ' 

Definition 2.1. For each n 6 N, we define a diffusive matrix process (Mn(t))t^ such 
that M^(0) := Mq and fort^O 

dM£(t) = -jM?(t)dt + e?dH? + (1 - e?)dY t (2.2) 

where (H^) t ^ is a d x d symmetric (resp. Hermitian) as defined in (|2.ip whereas 

d 



^ = v^£xr(M) dB\ 

i=l ^ ' 
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with i.i.d Brownian motions {B\)t^o and where Xi([ n t]/ n ) *s the spectral projector as- 
sociated to the i-th eigenvalue Xi([nt]/n) of the matrix Mn([nt]/n) if the eigenvalues are 
numbered as X\([nt\/n) < \2([nt]/n) < ••• < Xd([nt]/n) (we shall see that the above is 
possible as the eigenvalues are almost surely distinct at the given times {k/n,k 6 N}J. 

As for all t, the matrix Mn(t) is in the space we know that it can be decomposed 

as 

Mi(t) = 0^t)A^t)0^t)* 

where A^(i) is the diagonal matrix whose diagonal is the vector of the ordered eigenvalues 
of Af£(t) and where o£(t) is in the space for all t E R + . We also introduce a matrix 
O^(0) to be the initial orthogonal matrix (resp. unitary if j3 = 2) such that Mq (i) = 
O^(0)A O /3 (0)* where A := diag(Ai, X d ). 

The evolution of the eigenvalues of Mn(t) during the time interval [k/n;(k + l)/n] 
is given by independent Brownian motions if = and by Dyson Brownian motions if 

The eigenvectors of Mf[(t) do not evolve on intervals [k/n; (k + l)/n] such that = 
and evolve with the classical diffusion of the eigenvectors of Dyson Brownian motions if 
ejj = 1 (see [2] for a review on Dyson Brownian motion). 

Our main theorems describe the asymptotic properties of the ordered eigenvalues of 
the matrix Mn(t) denoted in the following as 

(A?(t) < m < ••• <A5(t)) (2.3) 

and also those of the matrix On(t) defined above, as n goes to infinity. 

Let (bl)t > Oj i £ {!)••• i d} be a family of independent Brownian motions on R. Recall 
that Cepa and Lepingle showed in [2] the uniqueness and existence of the strong solution 
to the stochastic differential system 

d\i(t) = -j\i(t)dt + V2db\ + /3p V 1 dt (2.4) 

starting from A(0) = (Ai ^ A2 ^ ■ ■ ■ ^ Ad) and such that for all t ^ 

Ai(t) < A 2 (t) < ••• < X d (t) a.s. (2.5) 
For the scaling limit of the ordered eigenvalues, we shall prove that 

Theorem 2.2. Let Mq be a symmetric (resp. Hermitian) matrix if (3 = 1 (resp. (3 = 2) 
with distinct eigenvalues X\ < A2 < • • • < A^ and {Mn{t))t>o be the matrix process defined 
in Definition \2.1\ Let A"(t) ^ ... ^ A^(i) be the ordered eigenvalues of the matrix Mn(t). 
Let also (Ai(i), . . . , Xd(t))t ^ ^ e the unique strong solution of (j2.4j) with initial conditions 
in t = given by (Ai , A2, . . . , A^) . 

Then, for any T < 00, f/ie process (A™(i), . . . , A^(t)) tg [ 0i T] converges in law as n goes 
to infinity towards the process (Ai(i), . . . , Ad(t)) tg [ 0i r] *n ^ e space 0/ continuous functions 
C([0,T],M. d ) embedded with the uniform topology. 
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In the case where j3p ^ 1, the eigenvalues almost never collide and we will see (see 
section 16. lh in this case that it is easy to construct a coupling of A and A n so that A n 
almost surely converges towards A. 

We shall also describe the scaling limit of the matrix On(t) (the columns of On(t) are 
the eigenvectors of Mn{t)) when n tends to infinity, at least until the first collision time for 
the eigenvalues, i.e. until the time T\ defined as T\ := inf{i ^ : 3i G {2, . . . , d}, Xi(t) = 
Vi(t)}. 

Let wf?j(t),l ^ i < j ^ d be a family of real or complex (whether /3 = 1 or 2) 

standard Brownian motions (i.e. w^(t) = B}-{t) + \^-T(/3 — l)Bfj(t) where the B^Bfj 
are standard Brownian motions on R), independent of the family of Brownian motions 
{b\)t^o,i £ {1, ...,d}. For i < j, set in addition Wj^t) := w^(t) and define the skew 
Hermitian matrix (i.e. such that R 13 = —(R 13 )*) by setting for i ^ j, 

a dwf- (t) a 

Then, with Aj(i),0 ^ t ^ T\,i £ {1, . . . ,d} being the solution of (|2.4p until its first 
collision time, there exists a unique strong solution (O^(t))o ^t^Ti to the stochastic dif- 
ferential equation 

dOP(t) = ^0?(t)dRP{t) - ^{t)d((R fi )*,R p ) t (2.6) 

This solution exists and is unique since it is a linear equation in and R^ is a well 
defined martingale at least until time T\. It can be shown as in [21 Lemma 4.3.4] that 
0^{t) is indeed an orthogonal (resp. unitary if f3 = 2) matrix for all t G [0;Ti]. 

We mention at this point that the matrix On(t) is not uniquely defined, even when 
we impose the diagonal matrix to have a non-decreasing diagonal A™(i) ^ ... ^ A n (t). 
Indeed, the matrix On(t) can be replaced, for example, by —On(t) (other possible matrices 
exist). The following proposition overcomes this difficulty. 

Define T n (l) to be the first collision time of the process (A™(i), . . . , A^(i)). 

Proposition 2.3. There exists a continuous process (On(t))o < t < Ti in with a uniquely 
defined law and such that for each t £ [0;T n (l)], we have 

0^t)A^t)0^tT l =Mi(t), 
where A^(i) is the diagonal matrix of the ordered (as in (j2.3j) ) eigenvalues of Mn(t). 

Proposition 12.31 is proved in Section [JJ We are now ready to state our main result for 
the convergence in law of the matrix On(t). 

Theorem 2.4. Let rj and T be positive real numbers. Then, conditionally on the sigma- 
algebra generated by (A"(s), . . . , A^(s)), 

^ s ^ T\ A T, the matrix process (On(t)) ^ t ^ m-^AT introduced in Proposition \2.3\ 
converges in law in the space of continuous functions C([O;T],0^) towards the unique 
solution of the stochastic differential equation (|2.6p . 



5 



Theorem 12.41 gives a convergence result as n goes to infinity for the eigenvectors of 
the matrix process (Mn(t)) but only until the first collision time T\. If p/3 ^ 1, the result 
is complete as one can show (see [2] and section [6J]) that the process (\%(t), . . . , \d(t)) 
is a non colliding process (i.e. almost surely T\ = oo). However, if p/3 < 1, it would be 
interesting to have a convergence on all compact sets [0; T] even after collisions occurred. 
Our next results describe the behavior of the columns of the matrix 0^(t) denoted as 
(01 (t), . . .,(f>d(t)) when t -> Ti with f < Ti. 

We first need to describe the behavior of the eigenvalues (Ai(i), . . . , (t) ) in the left 
vicinity of T\ . 

Proposition 2.5. If p/3 < 1 then almost surely T\ < oo and f/iere exists a unique index 
i* G {2,... , d} suc/i i/iai Aj*(Ti) = Aj*_i(7i). While we have, for all t and almost 
surely, 



The first part of Proposition 12. 51 is proved in subsections 13.11 and 13.21 the last statement 
is proved in [7J Hence equality (j2.7|) implies the existence of diverging integrals in the SDE 
(12. 6ft . Because of this singularity, we will show 

Proposition 2.6. Conditionally on (Ai(i), . . . , A^(i)), ^ t ^ Ti, we /lave: 

1. For all j ^ i*,i* — 1, the eigenvector <fij(t) for the eigenvalue Xj(t) converges almost 
surely to a vector denoted <j)j as t grows to T\. The family {4>j,j ^ — 1} is an 
orthonormal family of W 1 (respectively C d ) if f3 = 1 (resp. (3 = 2). We denote by 
V the corresponding generated subspace and by W its two dimensional orthogonal 
complementary in M. d (resp. C d ). 

2. The family {4>i* (t), <^*_i(i)} converges weakly to the uniform law on the orthonormal 
basis ofW as t grows to T\ . 

The paper is organized as follows. In Section [3l we review and establish some new 
properties for the limiting eigenvalues process (Ai(t), . . . , A^(i)) defined in 12.41 that will be 
useful later in our proof of Theorems 12.21 and 12.41 We also introduce, in subsection 13.41 
a process with fewer collisions that approximates the limiting eigenvalue process. In fact 
this gives a new construction of the limiting eigenvalues process already constructed in 
[1], perhaps simpler and more intuitive using only standard Ito's calculus. We give some 
useful estimates on the processes of eigenvalues and matrix entries of M„ in Section 01 In 
Section [5l we prove the almost sure convergence of the process (A™, . . . , AJJ) to the limiting 
eigenvalues process (Ai, . . . , A^) until the first hitting time of two particles with a coupling 
argument. In Section [6l we finish the proof of Theorem 12. 21 by approximating in the same 
way the process (A", . . . , Aj) with the same idea of separating the particles which collide 
by a distance 8 > 0. At this point, it suffices to apply that the result of Section [5] to show 
that the two approximating processes are close in the large n limit. In Section [7J we prove 
Theorem 12.41 the last statement of Proposition 12.51 and Propositions 12.31 and 12.61 




the following divergence occurs almost surely 




(2.7) 
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3 Properties of the limiting eigenvalues process 



In this section we shall study the unique strong solution of (12. 4ft introduced by Cepa and 
Lepingle in [JJ. We first derive some boundedness and smoothness properties. In view of 
proving the convergence of A n towards this process, and in particular to deal with possible 
collisions, we construct it for pj3 < 1 as the limit of a process which is defined similarly 
except when two particles hit, when we separate them by a (small) positive distance, see 
Definition 13.61 

3.1 Regularity properties of the limiting process 

Lemma 3.1. Let A = (Ai ^ A2 ^ • • • ^ A^). Then there exists a unique strong solution 
of (|'2.4|) . Moreover, it satisfies 

• For all T < 00, there exists a, Mq > finite so that for M ^ Mq 



Furthermore, there exists a, Mq > finite so that for M ^ Mq and % 7^ j, we have 



Proof. The existence and unicity of the strong solution is [H Proposition 3.2]. 

For the first point, we choose a twice continuously differentiable symmetric function 
(j), increasing on M + , which approximates smoothly \x\ in the neighborhood of the origin 
so that (f)(0) = 0, x<j)'(x) > 0, \(j)'(x)\ < c and \<fi"(x)\ < c, whereas \<j>(x)\ > \x\ x |x| A 1 



P max sup \Xi(t)\ ^ M < e 



a(M-M ) 2 



(3.1) 



• For all T < 00, all i, j G {1, . . . , d}, i 7^ j 






For all t, we have Xi(t)(f>' (Aj(i)) ^ 0, and also 



d 



^(Aj(t)) _ 1 
Ai(t)-Aj(t) 2 



d 



^(A t (t))-^(A,(t)) < d(d-l) 
Ai(t)-Aj(t) " 2 



EE 



EE 




00 



i=i 5 yt 



1=1 jr^i 



We deduce from the above arguments that there exists C > such that 
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By usual martingales inequality, as cf>' is uniformly bounded we know that, see e.g. [21 
Corollary H.13], 



P 



sup | V / 0'(Ai(t))dfei(i) |> M 

0<t<T ~[ Jo 



< exp(-— ) 



and therefore using the fact that |<^(a;)| > |x| x \x\ A 1, we deduce the first point with 

M o = I Ei=i + C r and « = 1/2CT. 

For the second point, we first remark as in the proof of [U Lemma 3.5] that for all 
i < d 

n f T dt n ^ f T dt 



o I A d (i) - \i(t) | I - Xj(t) 



j<d 



dt 



o X^-Xjit) 



= X d (T) - X d {0) - V2b d T + 7 / A d (t)dt . 

JO 

so that the first point gives the claim fo j = d. We then continue recursively. □ 
3.2 Estimates on collisions 

To obtain regularity estimates on the process A, we need to control the probability that 
more than two particles are close together. We shall prove, building on an idea from Cepa 
and Lepingle [5], that 

Lemma 3.2. For r > 3 and I C {1, . . . , d} with \I\ = r, set 

s l= D(M*)- A i(*)) 2 - 

We let, for e > 0, 

rj := inf{t > : min S T t < e} 

\I\=r 

Then, for any T > and n > 0, for any r > 3 there exists e r > which only depends 
on {Sq, \I\ > 3} so that 

P«<r)<r?. 

Proof. The proof is done by induction over r and we start with the case r = d, 
I = {1, . . . ,d}. Then, S verifies the following SDE (see e.g. Theorem 1]): 

dS t = -2-yStdt + AVd^/~S t dp t + adt 

where /3j is a a standard brownian motion and a = 2d(d — 1)(2 +p/3d). The square root of 
Pt '■= \fS~t verifies the SDE 

dpt = -IPtdt + 2Vdd(3 t + (^-2d) — . 

2 pt 
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In particular, one can check that, if a = 2 — ^7 = 2 — (d — 1)(1 + pj3d/2) 



dp? 



-arypt 



dt + 2\fdap?- 1 dl3 t 



Thus, as a < for d > 3, for any e > 0, p? Ard is bounded so that J p^ ATd d/3 s is a martingale 
and therefore 



T 



E[p? Ard ] <P%~a 7 E[p? ATd ]dt 
Jo 

By GronwalFs lemma, since sup 4 E[p" Arti ] is finite, we deduce that 

E[/£ AT „] < pff(l - — )e"^ + ^. 
M ' s «7 07 

As a consequence, since a < 0, we have 

e «/2 p(r c* < r)< E[5 a/2j = E[p «^ ^ ^ (1 _ J_ )e -a 7 T + rff 



«7 



07 



We can take e small enough to obtain the claim for r = d. 

We next assume that we have proved the claim for u ^ r + 1 and choose e r +i so that 
the probability that the hitting time is smaller than T is smaller than r//2. We can choose I 
to be connected without loss of generality as the A* are ordered. We let R = mm{r/ , tJ^} 
when r/ is the first time where S 1 reaches e. Again following [5], we have 



/•TAR 

log S T AR = log5 / -27T + 4 v / 2 £ / 

r JO 



Aj(t) - A fc (t) 



/•TAB 

EE/ 



k,jel ' 

Aj-(t)- A fc (t) 



5/ 



dti 



1 



1 



j,fce/ Z£J 



si 



Xj^-X^t) A fc (t)-A,(t) 



n/9 rTAR j f 

-4r[(r-l)(^r + l)-2] jf JL 



(3.2) 



Note that M t = 4^2 £ fcjje/ JtT* Aj(a) gJ A * (a) db^ is a martingale with bracket A t = 16r |f . 
For r > 3, 4r[(r — l){rp/3/2 + 1) — 2] > 2pf3 > and therefore we deduce 



E[logS£ Ai J > log 5 7 - 2 7 T + 2/3 P E 



tar 



dt 
Sj 



+ E 



^EE 

j.fceJ 10 



TAR 



Aj(t)-A fc (t) , 



1 



L Ai(t)-A,(t) X k {t)-Xi{t) 



}dt 



For j, /c E I, we cut the last integral over times 



n 



3,k 



{t<T/\R:Y^ 



1 



Z£j 



1 1 

Xj(t) - Xt(t) X k (t) - Xi(t) - sj } 
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so that 

_ v f (A,(t)-A fc (t)) 2 v 1 _ dt_ 

fi-,k, si ^W*)-AK*))(A*(t)-Mt)) J - y 5/ 

This term will therefore be compensated by the third term in (|3.2p . For the remaining 
term, if I ^ / is such that min^g/ |Aj — Aj| < mm^j |A& — Aj| for all k ^ I then if t G £1? k 
and i* 6 / is so that minj g / |A/ — Aj| = |Aj — A$*|, we get 

d — r 1 

> 



(A/(t) - \i*{t)) 2 ~ Si 
and therefore on rl + \ > t, 

e r+l < Si + ^(Xjit) - Xi(t)) 2 < Si + 2r(A*(t) - Mit)) 2 + 2Sl < (3 + 2r(d - r))S/ 

As a consequence, we have the bound for all j, k £ I, all t E £l c - t < R, 
\j(t)-\ k (t) 



SI 



^ - 1/JSI ^ - y/3 + 2r(d-r)/y/. 



which entails the existence of a finite constant c so that 



r Xj{ t)-X k (t) 1 1 ](] , 

kimh, $ ^)-Ht) x k (t)-Ht) 1 

c \ - x - l' T dt 

^ ~v^hh'° i A i(*)-A«(t) i ■ 

Using Lemma 13.11 we hence conclude that there exists a universal finite constant c' de- 
pending only on T so that 



E[logSf AR ] > log^-2 7 T 



\/ e r+l 

On the other hand, we have 

E[logS£ AjR ] < P(r/ < T) log(e) + E[ sup log 5/] 

0<t<T 



where the last term is bounded above by (13. ID . We deduce that 



|log(e)| ^rpi|log(e)| | log(e) | | log(e) | 

We finally choose e small enough so that the right hand side is smaller than rj/2 to 
conclude. □ 
We next show that not only collisions of three particles are rare but also two collisions 
of different particles rarely happen around the same time. 
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Lemma 3.3. For all i,j such that i + 1 < j, set 

r £ V = inf{t > : (A*(i) - A,_ x (t)) 2 + (A,-(t) - A^t)) 2 < e'}. 
Then, for any T > and r/ > 0, t/iere exists e' such that 



Proof. Using Ito's formula, it is easy to see that 

d ((Ai - Ai_!) 2 + (A j - Aj_i) 2 ) = 8(1 + p/3)di 
- 2 7 [(A, - Vi) 2 + (Xj - A,_i) 2 ] dt 

+ 2V2 [(A, - Ai_i)(e»; - dVf 1 ) + (Aj - A i _i)(^' - d&T 1 ) 

(Aj — Aj_i) 2 t ^ (Aj — Aj_i) 2 



2p0 



E 



' . (Aj — Afc)(Aj_x — A fc ) 
1,1 



+ E 



. (Aj - A fe )(A 



A fc ) 



(it. 



Set X 4 := (Xi(t) — Aj_i(t)) 2 + (Xj(t) — Aj_i(i)) 2 and note that the quadratic variation of 
(A, - Aj-iKrfjj - tifer 1 ) + (Aj - Aj-xXaVs - dbj' 1 ) 





is 2i. Thus there exists a standard Brownian motion B so that 
dX t = 8(1 + p0)<ft - 2jX t dt + A^ t dB t 

(Aj — Aj_i) 2 \ - (Aj — Aj_i) 2 



2p$ 



E 



+ E 



k ^ hi ( Ai ~ A fc)( A *-i ~ A fc) fc^ij ( A J ~ A fc)( A i-i _ A fe) 



(it. 



Note that, by the previous Lemma I3T21 we can choose e such that 

p[t £ 3 <tk|. 



(3.3) 



Moreover, for all t r| such that X t e/4, we have for all k ^ i — 

(A*-A fc )(Ai_i-A fc )(t) ^ §. 

o 

The same property holds for j. To finish the proof, we will use the fact that the sum in 
the last term is bounded for all t ^ r| such that X% ^ e/4. We thus need to introduce the 
process Y t defined by Y t = m.in(X t , f ). Let us set f{x) := min(x, e/4)~ pl3 . Note that / is 
a convex function M + — > M + and that the left-hand derivative of / is given by 

/i(x) = -p/3x-^- 1 l { ^ |} . 

Its second derivative in the sense of distributions is the positive measure 



f"(dx) =p(3 



p/3(pf3 + 1) 



r p/3+2 < f } 



dx . 
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Thus, by Ito-Tanaka formula, see e.g. [H Theorem 6.22], we have 



Jo 



{Xs < f } 



dX* 



+ o [PP 7 



e\-P/8-i 



+ 1) 

x pP+2 



L x t {X)dx 



where Lf(X) is the local time of X in x. By definition we have 

, P /3+2 Hx 8 <z}a{X,X) 







o XI 



and thus, we obtain 



1 



2//3 2 / y; 



P f3jY s - pl3 dt + 4Y S Pl3 ' dB & ^j 
((Ai - Aj_i)(s)) 2 



(3.4) 



+ E 



.A^M ((Ai " Afc)(s))((Aj - 1 " Afc)(s)) 
((A i -A i _i)( S )) 2 



Ix.< e /4d8 + -p0(-) Lt(X). 



,,_ ; _ 1)i ((Ai-A fe )( S ))((A i _ 1 -A fc )( s ))_ 

The definition of local time implies that, almost surely, Lf(X) ^ i. We thus deduce from 
(HHP that 



E 



TAt 1j ,At§ 



i 



e\-«8-i 



T + C 



E 



-p/3 



with C = (p(3j + 4p f3 (d — !)-)■ Gronwall's Lemma implies that 



E 



-p/3 



If e' < e/4, equation (13. 5h implies that 



e\-p0-i 



T exp(CT). 



T i t ^ T A r; 



< V^expCCT), 



(3.5) 

(3.6) 
□ 



Lemma 3.4. VKii/i i/ie same notations as in the previous Lemma \3.3[ we have almost 
surely 



(ST 

Taking e' small enough gives the result with (j3.3|) . 

As a direct consequence, we deduce the uniqueness of the i* of Proposition 12.51 



inf r n * 



+ 0O. 



(k,£):k+l<e 

In particular, this gives the unicity of the i* in Proposition \2_ 
Proof. It is enough to write that for all e > 

P ( fc mf £ r M < Tj < d 2 { max P (r M < T A r|) + P (r| < T) } 

and deduce from Lemmas 13.31 and 13.21 that the right hand side is as small as wished when 
e goes to zero. 

□ 



12 



3.3 Smoothness properties of the limiting process 
Lemma 3.5. We have the following smoothness properties: 

• For all T < co and e > 0, there exists C,c',c finite positive constants so that for all 
5, rj positive real numbers so that n < die 1 A 5e) we have 



max sup |Aj(s) — A,(i)| > 5 

l<i<d S < t <( 3+J? ) AT 3 
o < i < T 



For all T < oo and e > 0, i/iere exists C, c' , c finite positive constants so that for all 
5,n positive real numbers so that n < die 1 A be) we have 



du 

max sup / 

#j s <t<(s+r,)Ari Js WW-Xj{u)\ 
€ t < T 



> 5 



?7 V 



ce 4 /rj 



Proof. Let us first fix s G [0,T] and set / = {i G {2, . . . , d} : |A;(s) - Aj-i(s)| ^ e/3} 
and note that on the event {s ^ r|}, the connected subsets of / contain at most one 
element. Let T E = inf{i > s : infj^/|Aj(t) — Aj_i(t)| ^ £/4}. The continuity of the Aj 
implies that T e is almost surely strictly positive. 

If i £ I U {/ - 1}, then we have, for t G [s; (s + rj) A r e 3 A T E ] 



|Ai(i)-AiO)| ^ 7 / |Ai(«)|dM + ^|6j-6i|+p/3 f £ 

./ s J S ■/■ 



du 



. |Ai(n) - Aj-(w)| 

J 7^ 



< 7 / |Ai(u)|cfa + V2|6J - + 4p/3(d - 1) 



i - s 



Using (|3.ip and [21 Corollary H.13], it is easy to deduce that there exists a constant c > 
such that for n < eS/ {8p(3(d — 1)) 



max sup |Ai(i) — Aj(s)| ^5 

te[a;(s+f7)Ar|AT e ] 



< cde 2, > . 



(3.7) 



Now, if i G /, with the same argument as for (|3.7p (the drift term in the SDE satisfied by 
Aj + Aj_i is also bounded), we can show that there exists a constant c > such that 



sup |(Aj + Xi-i)(t) - (Aj + Xi-i)(s)\ ^ 5 

te[s;(s+r?)Ar|AT e ] 



On the other hand, the process Xi(t) := (Aj — Aj_i)(i) verifies 



< ce 2, < . 



(3.8) 



dxfit) = 4(1 + p/3)dt - 7x?(t)cft + 2 Xl {t)(db\ - db^ 1 ) 
_ 2p/3 V (M*)-Aj-i(*)) 2 * 
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The denominator in the last term of the above r.h.s is bounded below on the interval 
t e Is; (s + rj) A r| A T e ] by 2pP(d - 2)±. Thus, using again (pLE} and [2, Corollary H.13], 
we can show that for 5 > crj/e, 



sup \Xi (t)- Xi (s)\ ^ V5 
t£[s;(s+7?)Ar|AT e ] 



< 



sup |xj (t) — x i (s)\ ^ 5 

t6[s;(s+r?)Ar|AT e ] 



<C ce 2l ? 



(3.9) 



where the first inequality is due to the fact that x% is non-negative. Using (|3.8[) and (|3.9 
gives for r\ < 5e/c 



max sup — Aj(s)| ^ £ 

i£/U{/-l} te[s;(s+r))ATfAT e ] 



Thus, with (j3.7j) . we deduce that for 77 < 5e/c 

\K{t) - \i(s)\ ^6 
In particular, there exists d > so that if e 2 > crj, 



^ 2ccfe c 2" . 



max sup 

* te[s;(s+»?)Ar|AT e ] 



< 2cde °^ . 



[T 6 < (s + r/) A r £ 3 ] < 



max sup |Aj(i) — Aj(s)| ^ 5e/12 
* s<t<(s+r))AT e AT§ 



< 



AcdT „i^a 



which is as small as wished provided rj is chosen small enough. This allows to remove the 
stopping time and get for any fixed s < T, and 5 > crj/e 



max sup \\(t) — Aj(s)| ^5 
1 s<t<(s+J?)Ar| 



<2cde- c&4 /^ + 2dce- c ' £i ^. 



The uniform estimate on s is obtained as usual by taking s in a grid with mesh rj/2 up to 
divise 5 by two and to multiply the probability by 2T jr\. Thus we find constant c, d, and 
C so that if 77 < c(e 2 A 5e) we have 



max sup |Aj(i) — \{s)\ ^5 

1 s<t<(s + ri)f\T% 



< 9Z. ( e -cS 4 /2 v + e -c'e 4 / V 
~ 7] V 



The second control is a direct consequence of the first as we can first consider the cas 
j = d to deduce that for i < d 



du 



< \\d(t) - \ d (s)\ + V2\b d (t) - b d (s)\ 



Xd(u) - Xi(u) 

where the right hand side is continuous. We then consider recursively the other indices. □ 
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3.4 Approximation by less colliding processes 

When p(3 ^ 1, it is well known [2, Lemma 4.3.3] that the process A has almost surely 
no collision. In this case, the singularity of the drift which defines the SDE is not really 
important as it is almost always avoided. In the case p/3 < 1, we know that collisions 
occur and in fact can occur as much as for a Bessel process with small parameter. The 
singularity of the drift becomes important, in particular when we will show the convergence 
in law of the process of the eigenvalues A ra towards A. To this end, we show that A can be 
approximated by a process which does not spend too much time in collisions. 
For 5 > 0, we define a new process (Xf(t)) t ^ o as follows. 

Definition 3.6. LetTi := inf{i ^ : 3i / j,Xi(t) = Xj(t)} and for allt < T x , setX\(t) := 
Xi(t). For t > T\, we define the process recursively by setting for all £ ^ 2,X\(Tf) := 
Xf(Tg—) + id and for t > Tf, the process Xf(t) is defined up to time Tf +l := inf{t > Tf : 
3i j,Xf(t) = Xj(t)} as the unique strong solution of the system 

dxf(t) = -jxfm + V2dbt+ P pj2 x s m d \ Sm • ( 3 - 10 ) 



The main result of this section is that 



Theorem 3.7. Construct the process X with the same Brownian motion b. Then, for any 
time T > 0, any £ G (0,p/3/4) 

limP ( sup max . |A*(t) - Xf (t)\ < 5* ) = 1 . 

The theorem is a direct consequence of the following lemma and proposition. 

Lemma 3.8. Let 5 > 0. Construct the process X with the same Brownian motion b than 
X s . There exists a constant c > such that, almost surely, for all £ G N 

max sup \X\(t) - Aj(t)| ^ c5£ . 

1 ^ * ^ d s£ t s£ T/ 

To finish the proof it is enough to show that Tf goes to infinity for £ <C 1/5. This is 
the content of the next proposition. 

Proposition 3.9. Let T < oo, < £ < p/3/4 and L = [l/S 1 ^]. Then the probability 
P [Tf ^ rl vanishes when 5 goes to zero. 

Proof of Lemma 13.81 We proceed by induction over t to show that, for each £, 



^(Xf-Xifit)] 
\i=i J 



1/2 



sup [ > XXi-XiY(t) ) ^cS£ 

< t < Tf 



with c = (X)f = i i 2 = d(d + l)(2d + l)/6) 
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• We treat the case £ = 1. By definition of the processes, X s = A on [0, T-f). At time 
t = Tf, the separation procedure implies that 

d d 

X>? - A 4 ) 2 (Tf ) = £((Af - A 4 )(Tf -) + i5) 2 = c 2 5 2 . 

i=l i=l 

The property is true for £ = 1. 

• Suppose it is true for £ For i £ [T/,T/ +1 ), since A" 5 and A are driven by the same 
Brownian motion, we get 

d d 

- A,(i)) 2 = -2 7 ^(Af(t) - Ai(t)) a (ft 

i=l i=l 

Observe that 
i=l 



Af(t)-A|(i) m-Xjit) 



1 



= i E E (Af (*) - **) - (A,,) - A,(*))) 2 {m _^ {m _ m) 
< 

as the (Ai)x<j<d an d the (A|)i<j<<j are ordered. Thus, 

sup E( A '(*) - A *(*)) 2 < E( A W) - A ^ T /)) 2 - ( 3 - 12 ) 

In addition, because of the separation procedure at time Tf +1 , we have 

/ d \ V2 / d X 1/2 

f^(Af-A,) 2 (T/ +1 )J =(E((\'-^)( 3 ti-) + «) J 

/ 1 XV2 

< (E^-^^+i-)) +^^^+i) C , 

where we used the induction hypothesis in the last line. The proof is thus complete. 

□ 

Proof of Proposition 13.91 In the case p(3 > 1, it is well known [21 p. 252] that T\ is 
almost surely infinite and therefore the proposition is trivial. We hence restrict ourselves 
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to p/3 < 1. Let rj > 0. Let us define the stopping times 



T^ d := in£{t ^ : min sj' S < e} , 

|/|=3 



:= inf{i ^ : min ((Af - XUY + (A - A _ 1 )^)(t) < s}, 
1 ^ %i0 < d 

where 5/' := Yli jeii-^i ~ A^) 2 (i)- Set also := r|' 5 A r| ,<5 . We know from Lemmas 
and 13.31 that we can choose e small enough so that 

IP [4e A rl < T] < t?. 

The number e being fixed, it is then straightforward to see from Lemma 13.81 that there 
exists 5q small enough so that for all 5 ^ 5q, we have 



Now, we have 



^ Yl 1 {T e s +1 -T e s > Si} < T ; T £ > T t 



We need to show that the second term goes to when S — > 0. Let {J-i}t>o be the filtration 
of the driving Brownian motion. We will prove in Lemma 13.121 there exists a constant 
c > such that, on the event {r^ ^ T[}, almost surely 



Tf > 5* | T T s\ > ccT p/3 + f . 



In the following, we suppose that 5 is small enough so that c5 3> 6 p P +2 £ and 

5-Z T - 5~ p ^ ^ - 5- pl3+2 ^. For such 5, we have 



Tf > <5«} 



^{Tf + i~- ! 

!=1 
L 



T/ +1 - T/ > 6« | T T , < - r^+«; r e 5 > 



£=1 



TL X - If >$\F 



^p[t/ +1 -T/^;t^t1 



where we used the Tchebychev inequality in the last line. Using Lemma 13.101 we get that 
there exists a constant C > such that 



-4? 



which goes to when (5 goes to 0. The proposition is proved. 



□ 
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Lemma 3.10. Let £ G (0; 2). T/ien i/iere exists a constant C > suc/i i/iai ; almost surely, 

a* ^ r/ +1 - r/ | j- T J ^ CT (1 - p/3)(1 - 2 ' 1?) . (3.13) 



on ; t s £ ^ T[ 



Proof. We know that there are no multiple collisions nor simultaneous collisions (be- 
cause of Lemmas 13.21 and I3.3j) and therefore we can denote by i the unique element such 
that Xf (Tf-) = Af_ 1 (T/-) and (Xf - Af_ 1 )(T/) = 6. We have by Ito's formula 



dt 

+ 2p/3 



(3.14) 



(A?-A?-i)(f) 



(Af-AtJW 



-T_ 1 (Af-Ai)W(Af_ 1 -Af)(i) 



dt. 



Let us define the Bessel like process {X t )t > o by Xq = 5 and for t ^ 0, 

Li-l v ~ » dt 



dX t = V2{db, 



(3.15) 



Using the comparison theorem for SDE [8j Proposition 2.18] (note that the drifts are 
smooth before Tf +l — Tf), we know that for all t G [0,T* +1 — Tf), we have almost surely 



(\!-\U)(T e 5 + t)^X t . 



(3.16) 



Let us define := mf{t : Xt = 0}. It is clear that almost surely Tf +1 — Tf sj Tj^. 



We thus have on ^ T£ 



<^t/ +1 -t/|.f : 



We finally conclude using a classical result for Bessel process, see e.g. [?, (13)]; the density 
with respect to the Lebesgue measure on 



of the law of the random variable Tj^ is 



ps(t) 



1 



1 



n^)t \2t 



5 2 \ 2 _£ 

e at 



Hence we deduce that for £ < 2 there exists a constant c > such that 

For time £ G [0;T], we define the random set 

J t := {i€{2,...,d}:|Af-Af_ 1 |(t)<Vi/3}. 



□ 



(3.17) 



Note that, on the event f2 := {r^ ^ T}, for each t ^ T, the set It contains at most one 
element. For each £ G {1, . . . , L}, and i G {1, . . . , d}, we define the stopping times 

if(\/I/3) := inf{t > T/ : min|Af - X^Ut) > , 

j 

tj(i, y/e/6) := inf{t > T/ : min|A* - Af_ x |(t) ^ v^/6} . 



If i denotes the unique index such that Xf(Tf 
s \S 



mm^|Af-A* |(T/) > 



Ai_i(r/-), note that if T/ < r e 5 then 
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Lemma 3.11. If Tf ?S r| and if i denotes the (unique) index such that Xf(Tg- 
X.f_i(Tg—) , then there exists a constant c > and So > such that for all 5 ^ 6q, we have 



i7 s 



(3.18) 



Proof. Note that i is the unique element of the set I T s defined by (|3.17p for which 
| Af - Af_ 1 |(T/) = 5. For a = 1 - p/3 and t € [T/;T/ +1 ), we have by Ito's formula 



- A^H*) = - 7 a(Xf - A^H^t 
+ a(Af - XUr-'it) V2(dbl - db\~ l ) -pp 



(3.19) 



1 ( A i - A fc)(*)( A i-l - A fc)(*) 



(it. 



For 4 € [T/,T*], we deduce that 



d(Af - Af_j)«(t) > a(Af - XUr-Ht) V2{db\ - db^ 1 ) - c'(Xf - Xl^dt 
where d = a r y + f3p(d — 2)36/e. Let be the first time after Tf so that Xf — Xf' 1 reaches 

,rpS,K 

k<5. Then, as £ e+1 {Xf - >4-i) a ~ 1 ( t ) V2{dbl - dUf 1 ) is a martingale, we find that 



E 



(Af - Af_!) a (4(Vi/3) A y/i/6) A 2^) | T T s, K > 6 a exp (-</ T) . (3.20) 



Before time *£(*, ^e/6), (Aj - A J 5 _ 1 )(t) can not cancel if j ^ i. Therefore we can choose k 
small enough so that the last inequality implies 



E 



(Af - \Un4 (VS/3) A t?(i, v/i/6)) l {t j ( ^ /3)A t 1( i,^/6) < T/ +1 } I > ^ 6XP (" C ' T ) 



which can be rewriten using the fact that \Xf — Xf_ 1 \(tf(y/e/3) A tf(i, y/e/6)) ^ as 
follows 

pk 5 (v^/3)Atl(i,Vi/6)^r/ +1 | Jul - ' ^ 



> S" [ — ) exp(-c'T). 



The lemma follows with c = (-^) a exp(— d T) 



□ 



Lemma 3.12. Let £, T > 0. There exists a constant c > and <5o > so that if 5 < 5q, 
on T/ ^ r e 5 A T, 

P «S - T/ I J" Ta ] > cJ 1 ^ . (3.21) 



Proof. We assume in the sequel that 6 < 1. The proof is based on Lemma 13.111 It 
implies 



a* ^ t/ +1 - r/ i j; 



(v^/3) A tf(i, v/i/6) ^ r/ +1 ; ^ < T/ +1 - T/ < 1 I T % 
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By Lemma |3.11[ we deduce that 



< t; 



e+i 



Tf <1\ Tr 



T 



> cS 



l-p/3 



4(x/i/3) A if(i, v^i/6) sC T/ +1 < T + 1; 5^ T* - T/ | T : 



But 



P [4(^/3) A i 6 e (i, y/i/6) ^ T/ +1 < T + 1; T/ +1 - 7? ^ ^ | J" T| 
^ P \tf (Vi/3) sC A T + 1; T/ +1 - if (x/i/3) ^ <S* | J^- 
+ P fci, x/i/6) < 4(Vi/3); v^i/6) - T/ ^ ^ | J" T/ 
Let us handle the first term of the previous right hand side 
P k 5 (v^/3) < T/ +1 A (T + 1); T/ +1 - ff (Vi/3) < <5« | ^ ( ^ /3) 



max sup |A|(s) - Aj(t|(Ve/3)) | ^ | ^ ( ^ /3) 

J t!(v^/3)<s<(tf(V£/3)+<5«)Aif(V£/12)A(T+l) 24 



Cexp(--^-) 

where we used Lemma 13.51 for the last line (actually the proof since we used the estimate 
for a fixed s). For the second term, the idea is similar 



max 



sup 



|AS00-Aj(ZJ?)|£ 



j ^ T/ < s< (T/+5f)Af|(i >N /i/6)A(T+l) 



12 1 ^ 



sC Cexp( 



by Lemma 13.51 As for all £ > 0, exp(— <C p/3 for small enough 5, the proof is 
complete. □ 



4 Properties of the eigenvalues of M% 

In this section, we will study the regularity and boudedness properties of the eigenvalues 

Of Mn- 

Definition 4.1. Let Mq be a symmetric (resp. Hermitian) matrix if j3 = 1 (resp. j3 = 2) 
with distinct eigenvalues \\ < A2 < • • • < and {Mn(t))t ^0 be the matrix process defined 
in Definition ^. 1\ For all t ^ 0, the ordered eigenvalues of the matrix Mn(t) will be denoted 
6y A?(t) < A?(t) < ... <A5(i). 

The following proposition characterizes the evolution of the process X n (t) until its first 
collision time. 
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Proposition 4.2. Let (A™(i), . . . , A^(i)) &e the process defined in Definition \4-l\ and 
set T n (l) := inf{i ^ : 3i / j,Xf(t) = A™(i)}. Then, almost surely, the process 
(A"(i), . . . ,X2(t)) verifies for every k £ N, i/ie following strict inequality 

\ n 1 {k/n)<\ n 2 {k/n)<...<\ n d {k/n). (4.1) 

In addition, there exist a sequence of Bernoulli random variables (ejJ)fcgN with meanp and a 
sequence of independent (standard) Brownian motions (b\)t > o> i € {1, . . . , d} also indepen- 
dent of the Bernoulli random variables (e^)fceN such that, the process (A™(i), . . . , X%(t))t ^ o 
is the re-ordering of the process (/z™ (t), . . . , ^{t))t > o defined for t ^ by 

dtf(t) = - 7 ^(t) dt + >/2d&; + p y, n (i f_ Km dt • ^ 

iw'i/i initial conditions in t = given by (/i"(0), . . . , ^(0)) = (Ai, . . . , A^). /n particular, 
up to time T n (l), the process A n verifies 

n 

d\nt) = -i\m ^ + V2db\ +/3Y X n (t) l X n {t) ^ . 

jjLj i \ > j v J 

Remark here that we use the property that e" = (e™) 2 . 
Proof. Let us show first that for each k € N such that k/n < T n (l), we have almost 
surely the strict inequality (14. ip . We will proceed by induction over k. Note that under 
our assumptions, it is true for k = 0. Suppose it is true at rank k and let us show it is 
then true at rank k + 1. From Definition 12.11 if the eigenvalues of Mn(k/n) are denoted as 
\™(k/n) < • • • < X2(k/n), then, depending on the value of the Bernoulli random variable 
e£, the dynamic for t £ [k/n; (k + l)/n] is 

• if e2 = 1, the process (A"(i), . . . , A^(£)) follows the Dyson Brownian motion with 
initial conditions (A"(fc/n), . . . , A^(A;/n)) (see [2] Theorem 4.3.2]); More precisely, 
we have for t 6 [fc/n; (A; + l)/n) 

dA?(t) = - 7 A?(t) (it + v^dW? + p 



where the (W^ )i ^ o> * 6 {1> • • • , d} are independent Brownian motions. In particular, 
this process is non-colliding in the sense that the A"(i) will almost surely remain 
strictly ordered for all t G (/c + l)/ n ) (see [21 Theorem 4.3.2]). Thus, we will 

almost surely have A™((& + l)/n) < • • • < A^((A; + l)/n). 

on the other hand, if e£ = 0, we need to define a new process (//"(i), . . . , n2(t)) of in- 
dependent Ornstein-Uhlenbeck processes with initial conditions (A™(/c/n), . . . , Aj(A;/ti)); 
More precisely, the evolution for t £ [k/n; (k + l)/ra] is given by 

^(i) = -jtf(t)dt + \fldB\ (4.3) 

where the Brownian motions B % are the ones of Definition 12.11 Note that, before time 
T n (l), the two processes A n and \x n coincide. In this case, the fj,f (t) can cross and the 
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ordering can be broken in the interval [k/n; (k + l)/n]. However, if crossing for the 
process /i n happen before time t = (fe+l)/n still we know that e' y ^ k+1 ^ n fif((k+l)/n) 
are almost surely distinct. The re-ordering of the thus always gives ((fe+l)/n) < 
■■■<X n d {{k + l)/n)a.s. 

The induction is complete and proves equality (|4,ip for all k G N. We deduce from the 
above arguments that for k such that k/n < T n (l), the evolution of X n (t) for t G [fc/n; (& + 
l)/n A r n (l)) is 



dA?(i) = - 7 A?(t) eft + v^dW* + (1 - e?)dfl|) + /3 



5^) 



A?(t) 



di. 



with initial conditions in i 
cess 6* for t ^ by fe£ := 



= fc/n given by (A™(/c/n), . . . , \^{k/n)). Let us define the pro- 
J Q (e"dWj + (1 — e™)dB l s ). Using the fact that the Brownian 
motions (W^)t > o> i € {1, . . . , d} are mutually independent and independent of the Brow- 
nian motions {B\)f^o,i G {1, ...,d} (also mutually independent), it is straightforward 
to check that the processes (b\) t ^Q,i G {l,...,d} are mutually independent Brownian 
motions. It is also easy to see that, for all s,t G [k/n; (k + l)/«], the random variables 
e£(W| — W*) + (1 — <?1){BI — B l s ) and e£ are independent. Therefore, we deduce that the 
brownian motions (b l t )t ^ o,i G {1, . . . , d} are independent of the sequence (ej^fceN- □ 
The following regularity properties will be useful later on. 

Lemma 4.3. Let T < oo. T/ien f/iere exist constants C,Ao,c,c',a > which depend only 
on T, d such that for all n G N, all Aq and all £ > 



max sup |M^(t)jj| 
l < i,j ^ d < t < T 



max sup |M^(t)y 

1 < ij < d < s,i ^ T, 
|t-s| < 5 



> £ 



< Cexp(-a,4 2 ) 



< 7 ex P — -) . 

CO 



(4.4) 
(4.5) 



Proof. Using Ito's formula, we can check that 



•*MP{t)-e"Mfo) 



fe* UdHl + (1 - e^Y.xfi— 



Let us set A n (s,t) := e ft M^{t) - e 7S M^(s). The entries of A n (s,.) are martingales 
with respect to the filtration of the Brownian motions conditionally to the Bernoulli 
random variables (eiJ)fceN (this is due to the independence between the Brownian mo- 
tions {B\) t > 0) (Hf \ij))t ^ 0) 1 ^ hj ^ d and the sequence of Bernoulli random variables 
(ffc)fceN- Using the fact that |xr([ ns ]/ n )ijl ^ 1 f° r a ^ h3> we can check that there exists a 
constant C(d, T) which does not depend on n such that for all n G N 



\(A n (s,-)ij,A n (s,-) k i)t\ ^ C(T,d)\t-s\ 
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Let A > 0, using [21 corollary H.13], we have 



max sup |(e 7t Mf (i))^| > A 



< d 2 max P 

l<i,j<d 

= d 2 max I 

1 < i,j sc d 



sup |(e^Mf (t) - M$)ij\ > A - max |M£ (i, j)l 



sup |A„(0,%| > A-mBx\M§(i t j)\ 



.0 ^ t ^ T 



'J 



^ d 2 exp 



(A - max,,, | M{](i,j) \f 
C(d,T)T 



(4.6) 



Similarly, for any given s £ [0, T], for e > 0, using [21 Corollary H.13], we have, for each 
entry ij and for every 5 > 0: 



max sup |(e 7 *M^(i) - e 7S M^(s));.,-| > e 
1 s: i,j < rf( e [ s _5 >s+ 5] 



^ d 2 exp 



and therefore there exists a positive constant c' so that 



max SUp | V o xv± n yu) — c iv± n yo ))l3 

K w < a o < s,t < t, 

It— »| sj 5 



> £ 



[2T/d]+l 

i=i 
,2T 



max sup 

i<v<d| t _M | <(5/2 



|(e^(t) - e^/ 2 M£(ity2))«| > e/2 



< d — exp ( -— ) . 



□ 

Lemma 4.4. Let T < oo. Then there exist constants C, Aq, c' , c" , a, eo > which depend 
only on T, d such that for all n 6 N, all A ^ and all e > 



max sup |A"(t)| > ^4 



max sup |A"(t) - A"(s)| > e 

1 ^ i ^ d o ^ s,t ^ T, 

t-s| < 5 



s£ C"exp(-aA 2 ), 



c" , e\ 



(4.7) 
(4.8) 
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Proof. This lemma is a consequence of Lemma 14.31 and the inequalities 
max |AS(t)-Aj(«)| < fa |A?(t) - K(s)\ 2 ) 



vi=l / 



£ |M£(%-M<?( S ) y | 2 (4.9) 
^,i=i / 

<d max \M^t)ij - MP{s)ij\ 

±<%,j<d 

where, for the second inequality, we used O lemma 2.1.19] and the fact that the A™ are 
ordered. □ 

5 Convergence of the law of the eigenvalues till the first 
hitting time 

Proposition 5.1. Take A(0) = (Ai < A2 < • • • < A<f). Construct pJ 1 , strong solution of 
(|4.2p . with the same Brownian motion than A, strong solution of (|2.4|) . both starting from 
A(0). A n equals [i n till T n (l). For all T > 0, we have the following almost sure convergence 

lim max sup |A"(t) — \{(t)\ = 0. 



n— >oo 1 <g i < d 



t < TAT n (l)Ar^ 



As a consequence, if we let T\ = inf{t > 0, 3i 7^ j, Xi(t) = Xj(t)}, we have almost surely 

Ti < liminf T n (l) . 

We point out that this convergence does not happen on a trivial interval since we have 
Remark 5.2. For any -q > 0, there exists t(t/) > so i/ioi 

lim P[T n (l) >r(r?)] > 1 - rj . 

n— >oo 

Proof of Remark \5.SX By the same arguments developed in (|4.9p . we find that 



sup max |A™0)e 7 * - Aj(0)| ^ e 



< P sup |tr((M n (t)e 7 * - M ) 2 )| / e 2 

,2 

2 . 



^ d exp( 



2C(d,T)T' 



But since also the A" are uniformly bounded with high probability, we can choose for any 
r] > the parameter T small enough so that 



max sup |A™(i) — Aj(0)| > min |A, — Aj + i|/3 

l^i^d^<j! 1 < i ^ d 



This implies that P(T n (l) < T) < rj. □ 
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Proof of Proposition [5TT1 Using Ito's formula, we can compute 

£ (XUt) - Ht)) 2 = -Try f WOO " ^(«)) 2 ds 
i=i Jo i=i 

+2/3 r^x^Eww-M-)) 



i=i 



Af(s)-A?(s) Aj(s) — Aj(s) 



(5.1) 



ds 



By the same argument as in (13. lip the second term in the right hand side is non positive. 
Thus using equations 15. 1| we find for t ^ T n {\) 

Y (W) - mf < 2/3 - p) e E fil I A A 'LI da := • 

i=l J ° i=l j^j 1(1 ' ^ ' 

We next prove that 



lim sup -Rn(i) = a.s. 

n ^°° < t < TAr| 



(5.2) 



Write R n (t) as = P n (i) + Q n (t) where 



P„(i) := / ^-p)YY 



A?(H/n)-Ai( a ) 



<?»(*) == / k-p)EE 

70 i=i 



Ai(s) - Aj(s) 
i=i 



A?(s)-A?([ns]/n) 



Ai(s) - Aj(s) 



(is . 



We first handle the convergence of Q n (t). Set = {sup| s - t |<i/ n max!<j<d |A™(i)— A™(s)| ^ n 1 / 2+<E }. 

t<T 

On the event fii, we have 



|g n WI<-- 1/2+e EE / 

Following (|4.9p . we know that 



ds 



i=i i 



P(n?) < ce" 
We thus deduce from Lemma 13. II that 



I Ai(s) - Xj{s) 



sup |Q„(t)| > -5 

t < T 



ds 



5|^o |A,(s)-A i( s) 



> 



< ce 



+ ce 



Hence, Borel Cantelli's Lemma insures the almost sure convergence of Q n to zero. We now 
turn to the convergence of P n (t). Let 77 > small and write 



Pn(t) 



d(d - 1) 



(e™ - p)ds + P n (t) 
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ft * Xf ( [ns] In) - A™( [ns\ In) 

The process J *(e™— p)ds is a martingale and by Azuma-Hoeffding inequality, for any S > 



max | y (e™ -p)ds| > <A < 2exp( 



<5 2 n, 
'~2' 



We now use the independence between the brownian motions (b\)o ^ % ^ y, i = 1, . . . , d and 
the Bernoulli random variables e^,k = 1, . . . , [nT]. Conditionally on the (b l t )o <^t ^T,i = 
1, . . . , d, the processes \i(t), i = 1, . . . , d are deterministic and the process P n is a martin- 
gale with respect to the filtration of the e£. We let 

By Lemma 13.51 and Lemma 14.41 the set 

fl = { sup l^l^n- 1 / 8 } 

kKnThrf 

has probability larger than 1 — e -cnl/16 . Moreover, by martingale property it is easy to see 
that for all A > 0, 

E[l ne AP " (fc/n) -^ A2E ^ 1{ ^/" )2 ] < 1. 

Taking A = n 1 / 16 , since on tt, -n l / m \Al\ + n l l % \Al\ 2 /2 < 0, Tchebychev's inequality 
yields 

[Tn] \ 

{\P n (k/n A r 3 )| > n" 1 / 16 ^ \A%\ +f)}nfi < e~* 

£=0 / 

As by Lemma [3. 11 Y^i=o l-^&l * s bounded by n 1 / 32 with probability greater than 1 — e ~ nl/16 
we conclude that 

P (\P n {k/n A r 3 )| > n" 1 / 32 ) < Ce~ n ^ . 

The uniform estimate is obtained easily by controlling the increments of P n in between 
the times k/n,k < [nT] by sup fc< [ nT j \A%\ which we have already bounded. 

□ 



6 Proof of Theorem 12.2 , 
6.1 Non colliding case p/3 ^ 1 

It is straightforward to deduce Theorem 12.21 when p/3 ^ 1. Indeed if /?p ^ 1 we know that 
there are no collisions for the limiting process and more precisely, see e.g [21 p. 252], 

P(r 2 <T) ^c(A )T/|loge| 
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with some finite constant c(Aq) which only depends on the spacings of the eigenvalues at 
the initial time. This implies in particular that 

lim lim P(T" T) = 
from which we easily deduce Theorem 12.21 from Proposition 15.11 
6.2 Colliding case p(5 < 1 

We now define the process (A™ ,<5 (i))j ^ o which will depend on the sequence (Tf)f^ defined 
in Definition 13.61 To unify notations, set Tf := T\ and Tf(l) := T n (l). 

Definition 6.1. For t < Tf, set X™ ,S (t) := A™(i). For time t > Tf, we define the process 
recursively by setting for each 1^1, X™' S (Tf) = X^' S (Tf—) + i5 for all i G {1, ... ,d} 
and for t > Tf, the process A"' is defined up to time Tf +1 by ordering the process 

Gin*),..., 

^t^T 6 which is defined for t ^ Tg as 

— f n 

d^\t) = -jtf\t) dt + V2dbl + (3^ nS,, 1 nS,, dt ■ ( 6 - X ) 

with initial conditions in t = Tf given by (A™ ,l5 (T/), . . . , X^ ,S (Tf)). 

Lemma 6.2. Let T < oo and 5 > 0. We have the following convergence in probability, for 
all £eN, 

lim max sup | Xf (t) - A"' 5 (t)| = . 
"-*x>l<i<<*0<t<T/AT 

In particular, for every £, if Tf is the first collision time for X n,s after Tf_ x , 

Tf AT < liminf T*(£) A T a.s. 

Proof Again, we prove this Lemma by induction over i. 

• We begin with the case i = 1. Proposition 15.11 yields that the random variable 
maxi <^d sup <g t sg T n (i)/\T ~~ (*)l = converges to in probability as by Lemma 
13.21 P(r| > T) goes to one as e vanishes. Since we have the almost sure inequality 
Tf ^ liminfT^(l), the continuity of the Aj, 1 ^ i ^ d, the regularity property of the 
A™ given by Lemma 14.41 Lemma 13.51 and Proposition 15. 1} we can check that since before 
Tf Xf = Xi and Xf = A™, if Tf(l) < Tf A T, 

max sup |Af(t)-A^(t)| (6.2) 

l $ i $ ^ t<T?AT 

< max sup {|A? (t) - A"(T n 5 (l))| + |A,(t) - A,,(T*(1))|} (6.3) 

1 ^ 1 < "T«(l) ^ t<T{AT 

+ \xnrf(i)) - HTfm 

goes to zero in probability, when n goes to infinity. 
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• Suppose the property is true for £ and let us show that it is then true for I + 1. 
By the same argument as in the proof of Proposition 15.14 we can show that, for all t £ 
[T/; T*(l + 1) A we have 

E (K' 5 - a? ) 2 (*) < E (V - A 2 ( T /) ( 6 - 4 ) 

i=i i=i 

% ( ' p) ft *f (•)-*}(•)*' 

The same proof as in Proposition 15.11 shows that, if rf' is the stopping time rf for the 
process 



lim sup / f (eg -p)YY ^$ ds = o.s. (6.5) 

" — te[T/;r4(m)AT/+iAr 3.VT/ Af( s )-Af( s ) 



Thus, because of (|6,4p . the following convergence in holds 

lim max sup |A" ,<5 (f) - A?(t)| = o.a. (6.6) 

n ^°° 1 te[T/;T«(€+l)AT/ +1 A7f] 

Because of (|6.6p . we have T/ +1 Ar e 3 ^ liminf n _> 00 T*(i + 1) At £ 3 . Since the probability that 
rf is larger than T goes to one as e vanishes, we can show as in (|6.2p (note that Lemma 
14.41 Lemma 13.51 and Proposition 15. II extend to {A"' ,Xf,t > T/}) that in probability, 

lim max sup |Af(£) — X™ ,S (t)\ = 0. 

n^oo 1 < i ^ d T S( i+ y <; t <; T 6 +i 

The property at rank £ + 1 is established. The Lemma is proved. □ 

Lemma 6.3. There exists a constant c > such that for all L £ N, we have the following 
almost sure estimate 



max sup \\"' S (t)-\ 1 }(t)\^5Ly/c~. 



Proof. Note that the estimate is striaghtforward on [0, T-f ]. We then proceed by induc- 
tion on the time intervals [T/,T/ +1 ] as in the proof of Lemma 13.81 until the first collision 
time 

t! := mf{t > if : 3i, A?(t) = A?_i(t) or \f(t) = . 

We next claim that, at a given time, almost surely the eigenvalues A n are different. 
Indeed, this is clear if the eigenvalues follows Brownian motion and even more when they 
follow Dyson Brownian motion. Moreover the probability that more than two eigenvalues 
collide at some time vanishes. Indeed, this can only happen if the eigenvalues follow the 
Brownian motion. But the probability that 3 Brownian motions collide vanishes and hence 
the result. 
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Hence, there are almost surely at most two eigenvalues which can collide. Hence, 
let be the unique integer in {1, ...,d} such that A"(ii) = A™_ 1 (ti) (respectively 

A" ,<5 (ti) = X]' S (h)) and let n = ([rati] + l)/n. Notice that, for t £ [[rati]/ra; ([rati] + l)/n), 

we necessarily have e™ = 0. Let /^"' and /if for i G {1, . . . , d} be the processes such that 
for t G [ti; T\] 

d$(t) = - lfJ ,f(t)dt + V2db\ 

with initial conditions at t = ti respectively given by n n > s {t\) = X n ' S (ti) and fi n (ti) = 
X n (ti). We know that the A™' 5 , respectively the A™, are just a re-ordering of the processes 
tf' S and /if 

By definition, for t G [ti; ri], we find that : 

(^ _ M «)( t ) = e-TC*-*!)^ _ M » )(tl ) . 
As a consequence, we deduce that 

i=i i=i 

Moreover, as the A's are ordered but the set of the values of the A's and the ^'s are the 
same, using for instance [2j lemma 2.1.19], we have that 

j=l 3=1 

Gathering the above inequalities, we have shown that 

sup £(Af - X]) 2 (t) < Y,(\f - X"f(Tl) . 

*6[0,n] j=1 j=1 

We can continue inductively until we reach the time Tf +1 to finish the proof. 



7 Asymptotic properties of the eigenvectors 

Recall that ,i < j are real (respectively complex) standard Brownian motions if /3 = 1 
(resp. (3 = 2) with quadratic variation fit and that we also set for i < j, Wj i := w^-. In 
addition we also defined the skew Hermitian matrix R@ = — (R 13 )* by setting for i < j, 

a dwf- (t) a 

Proof of Proposition 



29 



It is classical to check that the unique strong solution of the stochastic differential 
equation 

n 

dO?(t) = e?OP(t)dR?(t) - j-OP(t)d((RPy,RP) t , (7.1) 

with initial condition On(0) := O^(0) (defined at the end of Section [Q), is in the space 0^ 
for all time t (see e.g. [21 Lemma 4.3.4]) and is such that, with An(t) being the diagonal 
matrix of the ordered (as in (|2.3[) ) eigenvalues of M^(i), we have 

0g(t)Ag(t)0g(*r l =M£(t). 

The law of the continuous process On is uniquely determined as the unique strong solution 
of ([71]). □ 
One can thus define the eigenvectors of Mn(t), denoted as (/>f(t), so that they satisfy 
the stochastic differential system 

«rw - «? E ^%>^> " f E ww f A?(()) ^ w M 

where to^-,i < j is a family of i.i.d. Brownian motions (on R if /3 = 1, C if ft = 2), 
independent of the eigenvalues A", 1 ^ i ^ ci. 



Proof of Theorem W. 



This proof is classical and uses the theory of stability for stochastic differential equa- 
tions. 

For T) > fixed, we deduce from Proposition 15.11 and Lemma 13.21 that the process 
(A™(i), . . . , X2(t)) converges almost surely in the space of continuous functions C([0; (T\ — 
if) A T],R rf ) (respectively C d ) if /3 = 1 (resp. f} = 2) endowed with the uniform norm 
towards (Ai(i), . . . , Ad(i)) ^ t ^ (Ti-^at where the Aj's are the unique strong solutions of 
(|2.4p (with the same Brownian motions b l ) and where Ti is the first collision time of the 
Aj, 1 ^ i ^ d. In the sequel we will work conditionally to the (A", Aj)'s satisfying the above 
convergence. 

Define for i / j the processes w^- n by setting 



B.n 



it) = J\yw^s) . (7.3) 



Note that the quadratic variation of this continuous martingale converges almost surely 
towards /3pt so that by Rebolledo's theorem (w^ n , i < j) converges towards (^/pw^,i < j). 

Moreover, if is the first time at which two eigenvalues are at distance less than e, the 
drift coefficients being bounded, we see, with a proof similar to the proof of Proposition 
15-H that for i ^= j 

rtAT( £ n 

-ds 



L 



o (\ n -A?) 2 00 



converges towards p Jq ATi (Aj(s) — Xj(s))~ 2 ds uniformly almost surely. Since T x e converges 
towards T\ as e goes to zero, the convergence holds till {T\ — rf) AT for any i] > 0. 

Gathering the above arguments, the result follows from [7J Theorem 6.9, p. 578]. □ 
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We now turn to the analysis of the behavior of the columns (f>i(t) of the matrix 0@(t) 
when t— >T\ with t <T\. Those vectors 4>i(t) form an orthonormal basis of M. d (respectively 
C d ) if {3 = 1 (resp. /3 = 2) and it is easy to check that they verify the following stochastic 
differential system 



#i(t) 



E 



dt 



2 ^(\ t (t)-\,{t)y 



(7.4) 



In the following of this section, we will denote by i* the unique (because of Lemma 
I3.4j) index such that Aj*(Ti) = Aj*_i(Ti). 

The main issue we meet at this point in the presence of collisions (that will occur if 
pj3 < 1; see [3]) lies in the divergence of the integral [2.71 that we now prove. 



We now describe the behavior of the d — 2 vectors (j>j(t),j ^ i* 



1 just before the 



first collision time T\. 

Proof of the first statement of Proposition \2.6\ 

We will denote by 4>je(t) the £-th entry of the <i-dimensional vector 4>j{t). For ^ t < Ti, 
we have 



dw? k (t)<f> k (t) - 



Xj(t)-x k (ty 



E 



(\j - x k y 



;(j)j(t)dt . 



(7.5) 



We recall from section 13.21 that there are no multiple collisions nor two collisions at the 
same time for the system (Ai(i), A2(t), • • • , A<j(i))o ^ t <Ti verifying (|2.4p . and therefore we 
may assume without loss of generality that for j ^ — 1, every diffusions and drift 
terms of (|7.5p remains almost surely bounded for t S [0;Ti]. To prove the lemma, we just 
need to prove that almost surely 



lim sup \\<f>j(t) - <f>j(s)\\' A 



0. 



s<T 1 



,i ; s ^ t<Ti 



The drift terms appearing in (|7.5p are obvious to deal with since l/(Aj — A&)(t) is bounded 
in the vicinity of T\ and that |0,-^(t)| ^ 1 for all t <T\. For the diffusion terms, we have 
for every I E {1, . . . , d} and for every s £ [0; T\] the following estimate 



sup 

«S t<Ti 



E 



Mi 



' Aj(u) - Afc(«) 



dw jk {u)(j) k i{u)\ > r, 



^ exp(- 



2fip(d - 1)M (Ti - s) 



) 



where M = sup ig [ 0;Tl ] max^j ^j- 
result. 



Using the Borel-Cantelli Lemma, we deduce the 

□ 

For 5 > 0, we want to define a process (4>i(t), (/>2(t), . . . , (pdity^-s < t<Tj that will be a 
good approximation of the process (4>i(t), (fe(t), . . . , </>d(i))Ti-<5 < t<Ti on the time interval 
[Ti — 5; T\\. Hence for j ^ i* ,i* — 1, we set <j>j(t) = 4>j (the vectors do not depend of time). 
It remains to define the evolution for (<^>j*_i(i), <^>j* (t)) that will depend of time t. 

Let V be the (d — 2)-dimensional subspace spanned by the orthonormal family {4>j]j ^ 
i* , i* — 1} and W its orthogonal complement in Let us define the "diffusive orthonormal 
basis" in the space W that will describe the evolution of the two vectors (<pi*-i(t), (pi* (t)) 
on the interval [Ti — 5; Ti] (up to the initial conditions at time t = T\ — 5 we will explicit 
later) . 
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Lemma 7.1. Let 5 > and (u,v) an orthonormal basis of the two-dimensional subspace 
W . We consider the following stochastic differential system 



(it 

T(V -Ai._i) 2 (t) 
p/3 dt 

T(v-V_i) 2 (t)' 



0i* (*) 



(7.6) 



d4>i*-\(t) 



(A;* - Ai._i)(t) 



(V - Ai*_i)(t) 



with initial conditions (0i*_i(Ti — 5),(f)i*(Ti — 5)) = (u,v). 

This stochastic differential system has a unique strong solution defined on the interval 



[Ti — 5;T\) such that for each t G [2j. — <5;7i), {(/>.;*_i(t), 0i* (t)} is an orthonormal basis 
ofW. 



Proof. For all e > 0, the function t — > l/(Aj* — Aj*_i)(t) is bounded on the interval 
[Ti — 5; Tf] and therefore there is a unique strong solution to the stochastic differential 
system (|7.6p till the time T-f where |Aj* — Aj*_i| < e as it is driven by bounded linear drifts. 
As Tf grows to T\ the proof is complete. 

To show that for all t G [Ti — S;T\) the family {0j*_i(t), (t)} is an orthonormal 
basis of W, we proceed along the same line as in the proof of [21 Lemma 4.3.4]. □ 

In the following lemma, we show that we can choose a constant 5 > small enough and 
an initial condition (u, v) G W such that the processes (4>i(t), . . . , (?!>i(t)) te [T 1 _5 ; j' 1 ) defined 
by Lemma ITTTI is indeed a good approximation of the process (0i(t), . . . , <^d(0)te[Ti-<5;Ti)- 
The advantage of the process (4>i(t), . . . , (/)i(t)) ig [j- 1 _5 : 2 i 1 ) is that it is simpler to study in 
the vicinity of T\ (see Lemma 17.31 below) . 

Lemma 7.2. Let ij > and k, > 0. Then there exists an orthonormal basis (u, v) of W 
and 5 > small enough such that if we denote by (</>j*_i(t), 4>i* (t)) te \x 1 -S;T 1 ) the unique 
strong solution of the stochastic differential system (|7.6p with initial conditions given in 
t = T x -5 by (0i»_i(t o ), 4>i* (to)) = (u,v), we have 



As for i G - 1} and j {i*,i* - 1} the terms l/(Aj - Aj) 2 (t) have almost surely 

a finite integral with respect to Lebesgue measure on the interval [tg;Ti) (in fact those 
terms are almost surely bounded as the corresponding particles remain at finite distance), 
the quadratic variation of the last term is of order 5 and therefore is smaller than n/2 with 
probability greater that 1 — k for 5 small enough. 

1 Note that all the diverging terms in Tj cancel in this expression. 
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It remains to check that we can choose (u, v) an orthonormal basis of W and 5 > 
such that 

\\fa*(T 1 -5)-u\\ 2 2 + \\fa*_ 1 (T 1 -8)-v\\ 2 2 ^ v /2. (7.8) 

This is a straightforward: Indeed we can approximate the 4>j{T\ — 5) for j {i* ,i* — 1} by 
the fa because of the first point of Proposition 12.61 thus we can choose two vectors {u, v} 
in the two dimensional space W so that (17. 8p holds. This completes the proof. □ 



We now turn to the study of the couple (fa*-i(t), fa* (t)) for t £ \T\ — <5;Ti) and in 
particular when t — ¥ Ti, t < T\. A crucial point is equation 12.71 which we now prove. 
Ito's Formula gives for t <T\ 



ln(Ai- \-i){t) = (- 7 + 2p/3)t + [ V2 

Jo 

-ppI* £ 



db\ — db % s 



(Aj* - Aj._i)(s) 
(is /"* 2 (is 



(Ai* - Aj)(V-i - Aj)(s) 7 (A^ - Ai»_i) 2 (s) 



If we suppose that J Q Tl dt/(\i* — Aj*_i) 2 (t) < +oo and since T\ < T e 3 for some e > small 
enough, we obtain a contradiction letting t — >■ Tj: under this assumption, the right hand 
side tends to — oo whereas the left hand side is almost surely bounded in this limit. □ 
The next Lemma [7.31 shows that the orthonormal basis (</>j*_i(i), fa* (t)) of the subspace 
W is in fact uniformly distributed in the set of all orthonormal basis of W in the limit 
t->Ti,t< Ti. 

As W is two dimensional, up to a change basis, we can suppose that the two vectors 
fa*—l(t) arid fa*(t) are two dimensional (we just study the evolution of their coordinates 
in an orthonormal basis of W). Let us define the two by two matrix fat) whose first line 
is the vector fa*(t) and second line is the vector fa*_i(t): 

fat) :-- 

Lemma 7.3. The matrix fat) converges in law when t —¥ T\,t < T\ to the Haar probability 
measure on the orthogonal group (respectively unitary group if (3 = 2.) 

Proof. To simplify notations, we do the proof in the case f3 = 1. 
Set to := Ti — 5 and define for t £ [0; 5) the function 

" t0+t ds 




<p(t) ■-- 



t0 (A*. - 

and denote by tp~ l its functional inverse. We now proceed to a change of time by setting 
for t e [0; 6) 

= fa*{tp- l (t)), &._i(t) = fo-i((p- x (t)) . 

As (p~ l {t) — > +oo when t — > S, t < 6 (because of (|2.7| ). the two by two matrix ip(t) whose 
first line is 4>i*(t) and second line is tpi*(t): 



i>At) 
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is now denned for all t £ R + and verifies the following stochastic differential equation 



vB ~ 

dil>(t) = y/pAi/}(t)dB t -^f^(t)dt. (7.9) 

where -B is a standard Brownian motion on R and where A is the two by two matrix 
defined by 

' r 



A - v-i o 

Note in particular that A 2 = —I. 

It is clear that there is pathwise uniqueness in the stochastic differential equation (17.91) 
(it is linear in ip). Therefore to solve entirely this equation, we just need to exhibit one 
solution. Using Ito's Formula, one can check that the solution is 



iP(t) = exp(^AB t )i>(0) 

cos(^/pB t ) sax(^/pB t ) 
-sin(y/pB t ) cos(^/pB t ) 



Note that for all t G R+, the matrix ip(t) is indeed in the space of orthogonal matrices. 

But (cos(y/pBt) , sm(^/pBt)) converges in law as time goes to infinity towards the law 
of (9, evl — 6 2 ) with 6 uniformly distributed on [—1, 1] and e = ±1 with probability 1/2, 
from which the result follows. □ 

Lemmas 17.21 and 17.31 give the second statement of Proposition 
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